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ABSTRACT 

7— I ■ 

i— i ' We have determined the mass profiles of dark halos in 83 objects observed by ASCA. 

■ The point spread function of X-ray telescope was deconvoled by the Richardson-Lucy 

■ algorithm and the temperature profiles were calculated to obtain the mass profiles. The 
derived mass profiles are consistent with the NFW model in 0.01 — 1.0r viria i. We found 
a good correlation between the scale radius r s and the characteristic mass density 5 C 
, which indicates the self-similarity of dark halos. The spectrum index of primordial 

S^-T density fluctuation, P(k) oc k n , was determined from the slope of r s — 5 C relation. For 

-^200 = 10 12 — 10 15 Mq, our analysis gives n = —1.2 ± 0.3 with a confidence level of 
90%. The mass density of dark halos is a good indicator of the mean mass density of the 
universe at the time when the halos were assembled, z = zj. Assuming 5 C oc (l+z^) 3 , we 
have determined the epoch when each dark halo was assembled. Our analysis indicates 
^ . that the field elliptical galaxies and groups of galaxies formed approximately at 1 + Zf ~ 

15 and at I + zj ~ 7 — 10 respectively. 

Subject headings: cosmology: dark matter — cosmology: observations — galaxies: clus- 
ters: general — galaxies: formation — galaxies: halos — X-ray: galaxies 



1. Introduction 

Early on, the density fluctuation grows linearly as the universe expands, dp/ p oc (1 + z)^ 1 . 
Once a density enhancement of a spherical region approaches dp/ p ~ 1, the sphere begins to turn 
around from the Hubble flow and collapse rapidly to form a virialized halo. The over density of 
virialized halo reaches 178po(l + z c ) 3 , where z c is the collapse redshift and po is the critical mass 
density in the present (Kaiser 1986). As the universe expands, the virial radius expands gradually 
by steady mass accretion. N-body simulations show that the accreting dark matter forms an 
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extended envelope, while the density profile of the original halo remains unchanged. The mass 
profile of a dark halo is rearranged only by a major merger through violent relaxation. A halo 
formed by a major merger at z rn is characterized by the mass density of 178po(l + z m) S , and it is 
preserved until the next major merger takes place (Salvador-Sole, Solanes & Manrique 1998; Raig, 
Gonzalez-Casado & Salvador-Sole 1998; Henriksen & Widrow 1999). By this means, the density 
profile of a dark halo can be a good indicator of the redshift when the last merger took place. 

Navarro, Frenk & White (1996) (hereafter NFW) suggested that all of the dark halos have the 
self-similar mass profile described as 

_ Po<>c 
P{r) ~ (r/r s )(l + r/r s )2' 

where po,r s and S c are the critical mass density of universe, the scale radius and the charac- 
teristic density respectively. In this case, the <5 C would be a direct indicator of the formation epochs 
of dark halos. 

In this paper, we will present the mass profiles of dark halos over the mass range of 10 12 — 
1O 15 M0 using a sample of 83 objects observed with ASCA. This furnishes us with a possible way 
to determine the formation epochs of galaxies and clusters of galaxies. 

Throughout the paper, we assume fi = 1, A = and Hq = 75kms~ 1 Mpc~ l . 

2. SAMPLE AND ANALYSIS 

Our sample consists of rich clusters, poor clusters, groups of galaxies and elliptical galaxies. 
The spectrum analysis was performed by the standard method using the Xselect and the Xspeck 
software to determine the HI absorption, the metal abundance and the luminosity weighted tem- 
perature. The results of spectrum analysis are described in the forthcoming paper along with a full 
description of our sample objects (Akimoto et al. 2000). 

The point spread function (PSF) of the ASCA telescope is characterized by a sharp central 
peak (FWHM ~ I5arcsec) and a significant scattering tail extending to ~ 3arcmin. We have 
deconvolved this scattering tail by the Richardson-Lucy algorithm (Lucy 1974). The observed 
images of 3C273 were employed as the PSFs. The deconvolutions were performed in the three 
energy bands (0.5 — 1.5, 1.5 — 3.0 and 3.0 — WkeV) separately, since the PSF has a weak dependence 
on the photon energy. All of the observed images went through 100 iterations. 

We have calculated the X-ray brightness profiles from the deconvoled images, masking the 
bright sub-peaks and obvious sub-clusters. The brightness profiles obtained in the lowest energy 
band, kT = 0.5 — l.5keV, were compared to those from ROSAT PSPC in the energy band of 
kT = 0.5 — 2.0keV, and good agreements were found in these independent observations. 
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The temperature profiles are calculated from the brightness ratio of the X-ray profiles in the 
three energy bands, kT = 0.5 — 1.5, 1.5 — 3.0and3.0 — lOkeV . We assumed that there were no radial 
gradients in metal abundance and HI absorption. The typical temperature profiles of our sample 
are shown in Figure 1, together with a comparison with the previous measurements. 



3. CALCULATION OF MASS PROFILE 

The total mass of a dark halo can be calculated from the gas density and the temperature 
profiles of ICM, assuming the hydrostatic equilibrium. 

If the observed X-ray surface brightness profile is described by the (3— model, S{6) = 
Sq(1 + (6 /6> c ) 2 )~ 3/3+0 - 5 , where 9, 9 C and (3 are the angular radius, the core radius and the (3 parameter 
respectively, and the radial gradient of temperatures is relatively small as compared with that of 
the gas density, the total mass within radius r is given by 

(r) G/^mp l + (r/r c ) 2 ' V ' 

where T, k, fj,m p and r c are the gas temperature, the Boltzmann's constant, the mean molecular 
weight of the hot gas and the core radius respectively. 

We found that r c , (3 and kT were not constant in a single cluster, but changed slowly with 
radius. The angular profiles of two parameters, r c tg^ and (3(0), were calculated by fitting the (3— 
model to the local brightness profiles of 9 =t 0.30. Although the angular profiles, r c ^\, (3iq\ and 
kT^, are the luminosity weighted properties of the hot gas along the line of sight, these are nearly 
equal to the radial profiles, r c ( r ), /?( r ) and kT^, if the radial gradient of gas density (or S^) is 
much larger than those of r c , (3 and kT. In our sample, r c ^ , (3^ and kTrm change factor 2 at 
most, while the changes 2 or 3 orders of magnitude. We therefore assumed that the observed 
angular profiles represented the radial profiles of these parameters. The local gas density profile 
was therefore given by nr r ) = uq{\ + (?V?* c (r)) 2 )~ 3/3(r) ^ 2 anc ^ t ne total mass Mi r \ was calculated from 
the equation (2). The mass density p( r +o.5dr) = (^(r+dr) — M^) / {ATrr 2 dr) was obtained from M( r ) 
and M( r+rfr ), keeping dr = 26arcsec. 

We have constructed the models of gas halos bound by the NFW dark halos to evaluate the 
systematic errors in our method. The mass profiles were calculated from the model profiles, which 
correspond to typical galaxies, poor clusters and rich clusters, by applying the same method. The 
best-fit 5 C and r s are then calculated by fitting the NFW model to the mass profiles. Comparing 
the calculated 5 C and r s with the original values, we have confirmed that the systematic errors of 
our method were much smaller than the typical photon noise of our sample objects. 



In our analysis, the brightness profile within 9 < 0.5arcmin were excluded to avoid the con- 
tribution of the cooling flow component. The maximum radii were extending up to 10 — AOarcmin 
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depending on the photon number available. We define r2oo as the radius where the mean interior 
density becomes 200 times the critical mass density of the universe, and M200 as the total mass 
within r2oo- In this paper, we consider r2oo and M200 is the virial radius and the virial mass (Cole 
k, Lacey 1996) respectively. 

4. THE OBSERVED MASS PROFILES 

The NFW model was fitted to our mass profiles to determine 5 C , r s , r2oo and M200 • The 
mass profiles of several objects, which have similar characteristic radius, are co-added to improve 
the signal to noise ratio. The composite mass profiles, scaled by r s are shown in Figure 2 by 
the shaded lines. The widths of lines correspond to the ±1<t errors. The average scale radii are 
< r s >= 24,59,150,200 and 460/cpc in (a), (b), (c) and (d) respectively. The NFW model and 
the (5— model were fitted to these composite profiles. Since the composite mass profiles drop more 
quickly at larger radii, the NFW model gives better fit to the observations. The x 2 /d.o.f. of the 
NFW model are 3.0/4, 4.3/5, 5.5/5 and 3.6/5 in (a), (b), (c) and (d), while those of the (3- model 
are 2.0/4, 24/5, 23/5 and 14/5. To illustrate a similarity of mass profile, we have normalized the 
composite density profiles by /?( rs )/10 9 and illustrated in the upper portion of Figure2 by the solid 
lines. It is remarkable that these haloes have very similar mass profiles. 

The brightness profiles scaled by /cT 1 / 2 (l + z) 9 / 2 can be a good indicator of mass profiles. Pon- 
man, Cannon & Navarro (1999) illustrated the scaled brightness profiles of 25 clusters, normalized 
by the viral radius, and pointed out the systematic change of profiles with kT l / 2 {l + z) 9 / 2 . Our 
measurements, however, do not conflict with this result. Since the observed dark halos have differ- 
ent concentration parameters c = r2oo/ r « with M200 as shown in the upper panel of Figure 3, the 
virial radius is not an adequate scale to illustrate the similarity of mass profiles. If our results are 
scaled by r2oo> we see the same systematic trend as found by Ponman, Cannon & Navarro (1999). 

The correlation between the M200 and the luminosity weighted temperature kT is shown in the 
lower panel of Figure 3. Our measurements give kT oc M 2 ° ^ 9±ai0 with a confidence level of 90%. 
We have found a good correlation between 5 C and r s in our 83 objects as shown in Figure 4. Our 
result is essentially the same as the c — M v i r and 5 C — M v i r correlations found by Wu & Xue (2000), 
but have larger dynamic ranges and smaller errors. The upper panel indicates residuals from the 
best-fit model. The rms scattering of S c shown by the dashed lines ( ±45% ) is significantly larger 
than the la errors of measurements. Our results indicate that the mass profile of a dark halo is 
described by a single parameter, such as r s , 5 C or M200. 
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5. THE SPECTRUM INDEX OF PRIMORDIAL FLUCTUATION AND THE 

FORMATION EPOCH OF GALAXIES 

The amplitude of primordial density fluctuation, P(k) oc k n , is characterized by the spectrum 
index n. In the CDM universe, the radius and the mass density of dark halo are connected to n as 
p oc r -(3«-+9)/( n +5) (Peebles 1980). The slope of r s — 5 C relation shown in Figure 4 is therefore the 
direct indicator of the power spectrum. If the P(k) is described by a single power in the mass range 
of our sample M200 = 10 12 — 1O 15 M0, the r s — 5 C relation gives n = — 1.2 ± 0.3 with a confidence 
level of 90%. The error is calculated from the intrinsic scattering of the objects. When the objects 
of M 200 < 2 x 10 14 M o are excluded, we get n = —0.8 ± 0.7(90%), which is consistent with the 
result of Wu & Xue (2000), n = -0.7 ± 0.3(la). 

We assume that the characteristic density 5 C is preserved during the steady accretion phase as 
N-body simulations suggest. In that case, S c should be proportional to the mean mass density of 
the universe when the halo was assembled, 5 C = f(l + Zf) 3 = /(l + z f,j) 3 [(l + zj)/(l + 2 ftj)] 3 , where 
/, 1 + Zf and 1 + z bj are a proportional constant, the formation epoch and the observed redshift 
of object respectively. As 1 + zj > 1 + z a bj, S c /(1 + z D 6j) 3 = /[(I + + z oty)] 3 > /• If our 

sample includes the zero age clusters Zf = z Q bj, the <5 C /(1 + -z D 6j) 3 of these objects would be equal to 
/. As shown in Figure 5, we found the clear locus in the lower end of (1 + z bj) v.s. S c /(1 + z bj) 3 
plot, and five clusters on the locus (A370, A1758, CL0016+16, MS04516 and RXJ1347.5-1145) 

were designated to the zero age clusters. These zero age clusters give the proportional constant and 

1/3 

the formation epoch as / = 1640 ± 260 and 1 + Zf = (0.085 ± 0.005)<5 C • Applying this equation 
to individual objects, we get 1 + Zf = 13 — 17 for the formation epochs of 3 elliptical galaxies 
in our sample (NGC1399, NGC3923 and NGC4636). Our results are nearly independent of the 
cosmological parameters as illustrated in Wu & Xue (2000). 

We are grateful to the ASCA team for their efforts on the design and operation of ASCA 
hardware and software. 
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Fig. 1. — The measured temperature profiles of ICM. The shaded boxes are our results. The 
height of the boxes indicates the error of 90% confidence level. The crosses are the results from 
Fukazawa et al. (1996) (NGC5044), Kikuchi et al. (1999) (MKW3S) and Markevitch (1998) (A3571 
and A3266). The vertical lines correspond to the 90% errors. The diamonds are from White (2000). 
The heights of diamonds indicate the la errors. 
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Fig. 2. — The composite mass profiles scaled by r s . The width of shaded lines indicates the la 
errors. The average scaling factors of emission measure are /cT 1//2 (l + z) 9 / 2 = 0.88, 2.3, 3.6, 6.1 and 
9.0 in (a), (b), (c), (d) and (e) respectively. The solid lines are mass profiles scaled by r s and 
/?( rs )/10 9 . The dashed line is the best-fit NFW model. The arrows indicate the virial radii r2oo- 
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Fig. 3. — The upper and lower panels are the measured M200 — c and M200 — kT relations. The kT is 
the luminosity weighted temperature determined by the X-ray spectrum. The error bars correspond 
to the 90% errors. The best fit relations with 90% error are logo = (7.2 ±1.8) - (0.44±0.13)ZosfM 2 oo 
and logkT = (-6.4 ± 1.4) + (0.49 ± 0.10)/o 5 M 20 o • 
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Fig. 4. — The measured S c — r s relation and the best-fit model. The error bars correspond to the 90% 
errors. The residuals are shown in the upper panel. The formation epoch of dark halos is indicated 
in the right side. The best fit relation with 90% error is log5 c = (7.61 ± 0.29) - (1.44 ± O.U)logr s . 
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Fig. 5. — A finding method for the zero age clusters. The (1 + z) is the redshift of each object and 
the <5 C / (1 + z) 3 indicate the age of dark halo since the last assembly. The dashed line is the average 
of 5 clusters on the lower boundary of 5 C /(1 + z) 3 — (1 + z) distribution. 



